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Abstract 

Electron stars are fluids of charged fermions in Anti-de Sitter spacetime. They 
are candidate holographic duals for gauge theories at finite charge density and exhibit 
emergent Lifshitz scaling at low energies. This paper computes in detail the field theory 
Green's function G^{uj,k) of the gauge- invariant fermionic operators making up the 
star. The Green's function contains a large number of closely spaced Fermi surfaces, 
the volumes of which add up to the total charge density in accordance with the Luttinger 
count. Excitations of the Fermi surfaces are long lived for uj < k^'. Beyond w fc^ the 
fermionic quasiparticles dissipate strongly into the critical Lifshitz sector. Fermions near 
this critical dispersion relation give interesting contributions to the optical conductivity. 
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1 Context 

Known large gauge theories with holographic gravity duals [Ij typically have a 'micro- 
scopic' field content of gauge bosons, scalars and fermions. These theories often have global 
U (1) symmetries under which a subset of the scalars and fermions can be charged. The the- 
ory can then be placed at nonzero chemical potential /x (equivalently, charge density (J*)) 
for this U{\) symmetry. The dual gravitational description of the nonzero charge density is 
that there must be an electric flux reaching out to the boundary of the bulk spacetime. The 
classical dynamics of the gravity dual will then determine the IR nature of the spacetime 
subject to this UV boundary condition. This general setup is illustrated in figure [T] below. 

The framework of the previous paragraph has similarities with difficult open questions in 
the theory of metallic quantum criticality in 2+1 dimensions. There, one typically considers 
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Figure 1: A nonzero charge density is imposed as a UV boundary condition. The IR physics 
is to be determined by solving the bulk equations of motion. 

a finite density of itinerant fermions interacting with gapless bosonic degrees of freedom 
[21 |3]. The strong effects of forward scattering of the infinitely many excitations at the 
(putative) Fermi surface invalidate perturbative frameworks for studying the low energy 
physics of these theories [HOE], although see also [ZllHlEj. Taking the ^7(1) symmetry of 
the large gauge theories as a proxy for the density of fermions (this identification will 
be literal if there are no massless microscopic bosons charged under the symmetry) we can 
use the holographic correspondence as a probe of the IR physics of metallically quantum 
critical states of matter. 

The presence of a density of fermions is a defining characteristic of metallic criticality. 
Fermionic spectral densities are therefore natural observables from which to understand the 
low energy physics. In particular one is interested in the possible smearing of fermionic 
quasiparticle poles at a Fermi surface due to interactions with the critical bosonic degrees 
of freedom. One should bear in mind that the fermions that can be accessed holographically 
are not the microscopic fermions but rather gauge-invariant composite fermionic operators. 

The simplest answer to the question posed by figure [l] above, that of determining the 
IR geometry sourcing an electric flux towards infinity, is the extremal Reissner-Nordstrom- 
AdS black hole jlOj . Here there are no charge-carrying fields in the spacetime and the flux 
emanates from a black hole horizon. The analysis of fermion spectral functions in extremal 
black hole backgrounds, initiated by led to interesting results [12l [TSj [H] . The fermion 
spectral functions exhibit Fermi surface singularities. Furthermore, if the gauge-invariant 
fermion operator has a sufficiently low scaling dimension relative to its charge, then the pole 
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is broadened into a 'non- Fermi liquid' branch cut singularity. The origin of this phenomenon 
was traced to the behavior of the Dirac equation in the far IR AdS2 geometry of the extremal 
black hole |14j and can therefore be thought of as indeed being due to strong interactions 
between fermionic excitations and the quantum critical bosonic modes in the IR. 

The chemical potential breaks Lorentz invariance and therefore an emergent scaling 
symmetry in the IR will generically be characterised by a dynamical critical exponent 
z 7^ 1 such that under scaling: t — t- At and x — )• X^^^x. It is natural to think of such 
emergent scaling in nonzero density holography as a strong coupling version of the Landau 
damping of bosons weakly coupled to a Fermi surface [2j. The Reissner-Nordstrom-AdS 
black hole is then seen to be an extreme case of Landau damping in which z — )• oo. It is 
of interest to consider holographic duals that admit, presumably more generic, finite z IR 
scaling behavior and to characterize the fermionic spectral functions in such theories. Some 
general observations about holographic fermions in finite z IR geometries were made in [T^] . 
In particular, it was noted that in a certain regime a Fermi surface pole would not only not 
be broadened, but that the width should be exponentially small at low energies, making 
the fermion excitations even more stable than those of a Fermi liquid! We shall reproduce 
this result, embedded within a slightly broader structure. 

Various interesting realisations of finite z IR geometries due to a nonzero charge density 
in the UV have recently been found. The first were in superconducting phases |16|I17]. Here 
the electric field is fully sourced by a condensed scalar field in the bulk, allowing the AdS2 
black hole throat to collapse into a 'Lifshitz' IR geometry |18tll9j. We are however interested 
in non-symmetry broken phases in this paper. A fluid of charged fermions outside of a black 
hole also absolves the horizon of needing to source the electric field and thereby allows a 
near horizon Lifshitz geometry [20]. The resulting geometries were dubbed 'electron stars' in 
|21j and have many features in common with the superconducting phases just mentioned. 
However, the U{1) symmetry is not broken by the fermions as these build up a Fermi 
surface in the bulk rather than forming a coherent state with a definite phase. Finally, 
without introducing charged matter explicitly, the extremal near horizon AdS2 throat can 
be collapsed by allowing a dilaton coupling to the Maxwell field in such a way that the 
effective Maxwell coupling goes to zero at the horizon [ 22l I23j . Our computations will 
focus on the electron star geometries for concreteness, and also because there is a certain 
minimalism in considering the spectral function of the same fermion that backreacts on the 
geometry. Nonetheless it will be clear that our results only depend on qualitative features of 
the geometry and will apply also to the other realisations of finite z we have just mentioned. 
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Sections 2 and 3 recall the essential features of electron star backgrounds and reduce the 
computation of the boundary field theory fermion Green's function to a quantum mechanics 
scattering problem. The fluid description of electron stars is the WKB limit of the bulk 
Dirac equation, and therefore the rest of the paper unfolds as a WKB analysis of scattering 
in a specific potential. Many of our results follow qualitatively from an inspection of the 
effective Schrodinger potential, as we discuss in section 3. 

Sections 4 through 7 characterise different low energy (oj <^ fj,) regimes of the fermion 
Green's function. The Green's function displays a large number of closely spaced Fermi 
surfaces with very long lived (as anticipated by ^15j) quasiparticle excitations. There is a 
maximal Fermi momentum kp. While this extremal momentum dominates the quantum 
oscillations of the star [24J, we show here it has the smallest spectral weight of all the 
Fermi surfaces. We show that the sum of volumes of the 2+1 dimensional field theory 
Fermi surfaces adds up to the total charge density, implying that the Luttinger count holds 
for this (non-Fermi liquid) system. Excitations of the Fermi surfaces become unstable at 
small momenta k < w^/^, where they dissipate efficiently into the critical Lifshitz sector. 
We characterise carefully the poles in the Green's function close to the critical dispersion 
u ^ k^ and estimate that they give a dominant contribution to the optical conductivity of 
the system. 

Finally in section 8 we take a limit in which the entire star is contained within the 
Lifshitz region of the spacetime. In this limit and with z = 2 we can solve for the Green's 
function exactly, thereby giving explicit examples of various statements made throughout 
the paper. 

2 The electron star background 

This paper will consider the fermionic spectral functions of field theories with electron star 
gravity duals. As mentioned in the introduction, our results mainly depend on qualitative 
features of the geometry that are shared by other spacetimes. We will emphasize these 
aspects throughout. 

The background spacetime and Maxwell field take the form 



with {f,g,h} functions of the radial coordinate r. In these coordinates r — )• at the 
spacetime boundary. Throughout, our spacetime will be 3+1 dimensional and the dual 




A 



hdt 



(1) 
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gauge theory, 2+1 dimensional. The equations of motion for the metric functions follow 
from the Einstein-Maxwell-charged fluid Lagrangian density 

as described in [21] . Here p is the local pressure of a bulk fermion fluid in terms of a local 
chemical potential and a local entropy density. The fluid is a corse-grained description of 
free Dirac fermions with unit charge (without loss of generality) and mass m. The fluid 
description is valid when [21] (analogous in depth discussions are also found in \lb\ [26] ) 

~ ^ « 1 , (3) 

and the following dimensionless mass is order one 

m = -^m . (4) 

This second condition implies that L?"im? ~ e^L?' / k? ~ 1/e^ ^ 1. The Compton wavelength 
of the fermion is much smaller than the radius of curvature of the spacetime and therefore we 
should expect to be in a WKB limit. This WKB limit will greatly simplify our discussion of 
the fermion Green's function below. While convenient, some of the essential physics should 
not depend strongly on this limit. We hope to explore a 'quantum electron star' explicitly 
in the future. 

In the deep IR, r — )• oo, the functions {f,g,h} are found to tend towards a Lifshitz 
geometry 



P 1 9c>0 , ^OO /f-s 



with in particular 



2 z-1 



hio = —- (6) 
The value of the dynamical critical exponent z depends on the order one dimensionless 
parameters m and c'^L'^/k^. It is bounded from below by [21] 

1 — 

This is the finite z emergent IR criticality advertised in the introduction. 

Integrating outwards from the IR, the local pressure of the fermion fluid decreases mono- 
tonically until it vanishes at a radius r = r^. This is the boundary of the electron star. 
Outside the electron star, r < r^, the spacetime is that of Reissner-Nordstrom-AdS 

+ g = -^, h = c[fi-rQ). (8) 

The four constants here {c, fi, M, Q} are determined by matching at r = and in turn 
determine the thermodynamics of the dual field theory. 




5 



3 The WKB Dirac equation 

Our objective is to compute the boundary Green's function of the gauge-invariant fermion 
operators that are making up the electron star. The bulk fermion can be taken to have 
four components and to satisfy the Dirac equation. With a judicious choice of basis for 
the Gamma functions, the Dirac equation separates into two decoupled equations for two 
two-component spinors. These equations differ only by the momentum k ^ —k. Without 
loss of generality we can focus on one of the two-component spinors. The Dirac equation 
can then be written 



y/g^asdr*^ -m^ - iy/g^ka2^ + + -4t)o-i$ = . (9) 

Here $ is a two component spinor, the cjj are Pauli matrices, and we have taken the field 
theory spacetime dependence $ cx e~*'^*+*'^^. The field $ is rescaled relative to the actual 
Dirac field: <I> = {—gg'^'^)^'^). 

The limits ([s]) and Q force us to consider the WKB limit of the Dirac equation. The 
large parameter that appears will be 



eL 

— = 7 > 1. 



(10) 



In particular, we can write the spinor as 

^=(^\Ae'^^. (11) 



Furthermore, in order to obtain an interesting uj and k dependence, we should rescale 

a; = 7^1) , k = ^k . (12) 

Then to leading order in 7 S> 1 the Dirac equations become 

iS' (h + u)^ ~ „ 

^ + - j^-rk(l)-m = 0, 13 

iS' (h + oj) rk ^ ,^ 
-— + ^^^ + --m = 0. 14 

These equations are simple to solve. The solutions are 

S = ±J dr^^ {uj + /i)2 // - r2fc2 - m2 , (15) 

TO- T * a/ {'^ + h)"^ I f — r'^k'^ — rh? 

cj) = ^ — . . (16) 

{u; + h)/^-rk 
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The overall function A is easily determined by working to next order in the WKB expansion. 
The zero temperature 'horizon' is at r — )• oo. We must impose ingoing boundary conditions 
at the horizon. The ingoing solution is, for oj > 0, the upper of the above signs. However, 
we must be careful to discuss the matching of the solution across the turning points of the 
exponent S. 

In order to match across the turning points it is convenient to map the Dirac equation 
into a Schrodinger form. This will allow us to use the standard machinery from WKB 



approximations in quantum mechanics. We will thus rederive the expressions (15) and (16) 
from this perspective. Writing 

it is simple to obtain a second order equation for each component of the spinor. The 
expression is a little clumsy, but it simplifies in our limit 7 ^ 1 to 

<^'l = I'g (m' + r^P - (^i±^) d,, , (18) 

$2 = ^fm$i -- — $;]. (19) 

(Cj + h)/^-rk\ 1^9 J 



Thus we recover the leading order exponents (15) in the usual WKB fashion. We will 



proceed shortly to perform a WKB study of the 'zero energy' Schrodinger equation (18) in 
the following section, but some comments first. The first of the two equations above has 
corrections of order 7]^ These corrections give a contribution to the WKB wavefunction 
of the same order as the standard V~i WKB prefactor, with V the Schrodinger potential. 
That is to say, a correction of order 7*^ in the exponent. Such corrections can have an effect 
on the phases that appear in performing matchings at points where the potential vanishes. 
These phases are responsible for the '^vr' type shifts in the ground state energy. In this 
paper we will not be overly concerned with the precise energies of the low lying states, but 
rather in states with intermediate to large excitation number. For these states, the effects 



of corrections to (18) are negligible (although computable if necessary) and we shall not 
consider them further. Note however that the final term in the second equation above needs 
to be kept in the WKB limit as the derivative brings down an extra power of 7. 



A second comment concerning (19) is that it appears to become singular when the 
denominator (a) + h)/\/~f — rk vanishes. The higher order corrections to the potential in 



^For completeness, the first order correction to the right hand side of equation (18 1, the Schrodinger 
potential, is 7my^ ^"^^^^:?'_^^^-^l"^^^~'' . In order to obtain this expression it is necessary to rescale $1 to 
eliminate a first derivative term that appears in ( |18[ ) at this order and thereby maintain the Schrodinger 
form. 
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(18) in fact also diverge at such points. These divergences are an artifact of ehminating 
variables and are not inherent to the Dirac equation Q. The easiest way in practice to 
avoid this complication is to note that if we instead obtain a Schrodinger equation for <I>2 
and then obtain <I>i from ^2, the divergence now occurs at {uj + h)/^/f + rk = while the 



effective Schrodinger equation remains the same (18). Thus to avoid possible singularities, 
we can swap between solving for or <I>2 first if we approach a singular point. We are 
interested in real k and r only, so nontrivial monodromies should not arise. The upshot is 



that we can just focus on the Schrodinger equation (18). We see that only k appears in 
the effective Schrodinger equation and so in the following we will restrict to > without 
loss of generality. 

The nature of the WKB solutions is controlled by the turning points in the exponent. 
From the turning points alone we will obtain a good picture of the Green's function we are 
after. Thus, we need to know the sign of the potential appearing in the effective zero energy 



Schrodinger equation (18) 



V = g [ r'^k'^ + rh 



2 £.2 , ^2 + 



f 



(20) 



In discussing this quantity, it will be useful to introduce the following notions. The local 
Fermi momentum is 

r^kUr) = j-m\ (21) 



The boundary of the star r = is determined by |21] 

kF{rs) = 0. (22) 

Note that kp{r) is negative for r < r^, i.e. outside the star. All the interesting features we 
are about to explore occur inside the electron star. The 'extremal' local Fermi surface is 
determined by ^24] 

kp = kp{r-i,) = max/ci?(r) . (23) 

The extremal radius always exists and is unique with > r^, inside the star. It can be 
determined numerically [21]. The local Fermi momentum kp{r) goes to zero at the center 
of the star as well as at the star boundary. 

We are interested in the low frequency Green's function oo ji. Start by setting w = 0. 
At zero frequency 

V = gr'^(k^ - Ppir)) , (w = 0) . (24) 
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Thus we have, for oj = 0, 

k > kp =^ V > everywhere , (25) 
k < kp =^ V < for ri < r < r2 and V > elsewhere . (26) 

Here clearly ri^2 are defined as the two solutfons to k = kF{r). Thus we see that the zero 
frequency WKB Dirac solution is exponential wherever the momentum is greater than the 
local Fermi momentum and oscillating whenever the momentum is lower than the local 
Fermi momentum. The oscillations therefore occur when the local states with the given 
momentum are filled. We may expect some resonant behaviour at k = kp separating the 
two regimes. The two cases are illustrated in figure [2| 




horizon horizon 

Figure 2: Illustrative radial dependence of the Schrodinger potential when cj = 0. The thick 
green line denotes oscillatory regions. To the left is the horizon while the boundary is to 
the right at r — >• 0. 



Now consider low but nonzero frequencies < |ct;| ^ ^. This should only affect the zero 



frequency results (25) and (26) in the near horizon region. This is because the near horizon 
(/, /i — )• 0) and zero frequency limits do not commute. Away from the infinite redshift of 
the horizon, the effects of small nonzero co are uniformly small. 
In the near horizon Lifshitz region ([s]) 



^Lif. = gooik + 



rh^ — (/loo + r^(^y 



(27) 



In the outer limits of the near horizon region, where cjr^ <C 1, we recover the zero frequency 



conditions of (25) and (26). This is the regime of overlap of the near and far regions and 
is consistent with our previous remark that small w should not qualitatively change the 
turning points in the far region. Note that the relations ^ and ([7]) imply that in the near 
horizon region lim r'^k'i, = /i^ — > 0. This is the consistency condition that the interior 
of the star is in fact occupied by fermions. 
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At any nonzero frequency we see that, assuming z > 1, for any mass and momentum, 
V < in the far interior r — )• oo. This fact may require another change in sign of V at 
a radius in the near horizon region and leaves us with three possible patterns of signs 
shown in the following table. The possibilities are also illustrated in figure [3] below. 



Case 


r — )■ oo 


r-i > r > r2 


r2 > r > ri 


r — ^ 


/ : k> k*p 


y < 


y > 


II : k<k*p 


y < 


V >Q 


V <0 


V>0 


III : k< k*p 


y < 


V>0 



Table 1: Three possible patterns of signs at a small but nonzero frequency uj. 



For the the third case shown in table [T] to occur, it must be that the second intersection 
point r2 at a; = occurred within the near horizon region. This is because the effects of 
finite uj are negligible outside this region. Furthermore, it must be that the momentum k 



Case I 





Case III 



1/r 



horizon 



horizon 



horizon 



Figure 3: Illustrative radial dependence of the Schrodinger potential (with z > 1) m the 
three possible cases at small nonzero frequency. The thick green line denotes oscillatory 
regions. To the left is the horizon while the boundary is to the right at r — )• 0. 



is small, so that its effects in making V positive can be offset with a small frequency uj. In 



fact, checking for the existence of zeros of (27) we can determine precisely when the third 
case can occur 

Casein ^ uj>\k\' max ^""^ ^ ~ . (28) 

0<a;<oo 

It is straightforward to find the maximum analytically, and we shall do so in a later section. 

In the following sections we will obtain the boundary fermion Green's function in some 
detail. However, a qualitative understanding of the frequency and momentum dependence 
of the spectral density - the imaginary part of the Green's function - can be extracted 
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simply from the Schrodinger potentials of figure [3| This gives us a global overview of what 
to expect in the following. 

We can see that in case II there will be many (because we are in the WKB limit) almost 
bound states that are stable up to an exponentially small tunneling amplitude into the 
horizon. These correspond to many Fermi surface poles in the retarded Green's function. 
At small frequencies the poles are bounded above in momentum by A; = A;^ and below by 
the Lifshitz dispersion k ~ uj^^^. These are the boundaries with cases I and III respec- 
tively. The boundary regions themselves are particularly interesting and we will consider 
them carefully below. In cases I and III we should expect no sharp features in the spectral 
density as there are no almost stable states. Putting these statements together leads to 
the schematic figure |4j This plot is very similar to plots describing the fermion spectral 
density in a superconducting background [281 [29]. This is unsurprising because the back- 




Figure 4: Schematic dependence of the dual field theory fermion spectral density as a 
function of uj and k for z > 1. Red lines denote the dispersion of sharp Fermi surface poles. 
In case I there are no low energy excitations because all the available Fermi momenta are 
too far away. In case III there are no well defined excitations because the Fermi surface has 
become unstable towards dissipation into the critical Lifshitz sector. 



ground fiuid of charged fermions has many features in common with a background charged 
bosonic condensate. Furthermore, the backgrounds [171 121] considered in [281 ISS] also had 
an emergent {z = 1) scaling symmetry in the IR and therefore a 'light-cone' within which 
fermion excitations were efficiently dissipated. One special feature of the z = 1 case is 
that the lightcone exists for both positive and negative frequencies. In our figure we have 
a 'non-relativistic lightcone' that only exists at positive frequencies. We will discuss the 
physics underlying figure [4] in the following section. In a later section we will present a 
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non-schematic version of this figure. We will also be concerned with the extent to which 
the sharp Fermi surfaces can be 'smeared' into a Fermi ball [21] and with the role of the 
relative residues of Fermi surface poles. Some of the Fermi surfaces will have a spectral 
weight that is exponentially suppressed relative to others. 



4 The boundary fermion Green's function: generic cj, k 

The Green's function of the fermion operator in the quantum field theory dual to our bulk 
fermion is read off from the asymptotic behaviour of solutions to the bulk Dirac equation. 
Dropping the subscript 1 for the spinor component, recall that we are interested in solving 



the Schrodinger equation ( 18 ) 



$" + 7V$ = 0, (29) 



with the potential given by ( 20 ) . Given that 7 1 , the WKB approximation will be valid as 
long as the potential and its first two derivatives remain of order one. This follows from the 
fact that the WKB approximation requires that the derivatives of the potential are small 
compared to the potential. We can use dimensional analysis to see how the comparison 
should be made. Because V has dimension two, we require: 

|F'(r)| < 7|y(r)|i , and \V" {r)\ <^ -f'^\V{r)\'^ . (30) 

These inequalities are obviously satisfied if V is order one with 7^1. Whenever V is 
smaller than order one we should consider a matching strategy over the problematic region. 
We will come to these cases later. 

In the outer region, defined to be the region from the boundary to the first turning 



point, we have V > 0. The two independent solutions of (29) are therefore schematically 



(7V)-3exp|±7^ dr\/y| . (31) 

Here e ^ 1 is a UV cutoff. Outside the star we must use the Reissner-Nordstrom-AdS 
expressions (|8]). It is then possible to show that to leading order near the boundary (recall 
mL > 1) 

<P = a+{uj,k)r"'^ + a-{u},ky-"'^ . (32) 

To obtain this precise result, including the leading correction to the powers itmL, one needs 
to include the leading correction, in footnote [T] above, to the Schrodinger potential as well 



as the V * prefactor in the WKB expression (|31|). This is the only place in this paper 



where we will refer to this footnote and indeed we do not even strictly need this subleading 
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correction here. Imposing ingoing boundary conditions at the horizon and matching across 
the turning points, we wih obtain some combination of these two modes in the outer region. 



Being careful to solve both (18) and (19) consistently and keeping the leading terms near 
the boundary r — )• 



Pouter = vl(w, A;)r— ^ + B{u;, A;)r"^ . (33) 
The retarded Green's function of the boundary operator dual to ^ is then (e.g. [H]) 



Because we will be focussing on the single Schodinger equation (29) for <I> it is convenient to 



re-express the Green's function in terms of the asymptotic behaviour (32) of this function. 



Thus, using (19), we obtain 



yRf, . _ A + t^/c — /c a+(a;. A;) 



Here \i and c are the boundary chemical potential and speed of light as appearing in ([8| . /t 
has slightly unconventional units due to compatibility with the notation in previous works. 



We now proceed to obtain the Green's function (34). Let us first consider the case in 
which V vanishes linearly in radius at all turning points. This is the generic situation in 
which we are not at a boundary between two of the three regions above. The matching 
computation is standard and we have relegated the details to appendix [A} The result for 
the retarded Green's function at leading order in the WKB expansion is, in cases I and III, 

G^{uj, k)oc^ lim r'^""^ exp | -2-f drVV^ , (cases III and I) . (36) 

This expression is of course precisely what we should expect. The Green's function is pure 
imaginary and given by the probability of the particle tunneling through to the interior 
oscillatory region and into the horizon. This probability is exponentially small. 

In case II the fermion propagates through an intermediate oscillatory region. The 
Green's function is found to be 

G^(a;, fe) oc ^ lim r'^"^^ exp | -27 ^ ' drVV^ , (case II) (37) 
cosh(X + iY) + sinh(X - iY) 



where 



cosh(X - iY) - sinh(X + iY) 



(38) 



where here 



X = -f drVV + log2, Y = -f dr^-V^-. (39) 
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The Green's function in case II has poles when ie~'^^ = tanY. Note that both X and 
Y are large positive numbers when u is real. Including the exponentially small imaginary 
part of the pole to leading order, the poles are at 

y = 7rn + ie-^^ . (40) 

Here n is an integer. This formula has an immediate interpretation. For a given momentum 
k < kp it is to be solved for the 'Bohr-Sommerfeld' frequencies These frequencies 
are closely spaced and have a real part much greater than their imaginary part. These 
resonant frequencies correspond to exciting a fermion that is making up the electron star 
background. The small imaginary part gives the probability for the excitation to tunnel 
through the barrier and fall into the event horizon. From a field theory perspective this 
is the probability of dissipating into the IR quantum critical modes. In the background 
electron star itself, this process is prevented through the electrostatic repulsion between 
electrons as well as the degeneracy pressure. 

Let us emphasize that there are two distinct exponential suppressions here. The first is 
common to all cases and is due to the fact that the Green's function is a boundary Green's 
function, that the potential is positive at the boundary, and that therefore the fermion 
needs to tunnel into the spacetime. This suppression is an overall factor, yet it will play 
an important role shortly as it means that some quasiparticle poles have a residue that is 
exponentially smaller than others, depending on how far the fermion has to tunnel. The 



second exponential suppression is the imaginary part of the poles appearing in ( 40 ) . This is 
the statement that the bound states in the star are very stable as they are separated from 
the horizon by a potential barrier. 

The next step is to translate the above results into explicit statements about the fre- 
quency and momentum dependence of the retarded Green's function. We are interested in 
low frequencies but arbitrary momentum. In this regime we have noted that there are three 
distinct cases 

k > kp Case I <-)• high momentum , 

k*p>k>Oj'/' o Case II o intermediate momentum , (41) 
^i/z > ^ ^ Case III -f-)- low momentum . 

See also figure [4] above. 

Case I is familiar in essence from Fermi liquids. At zero temperature and at momenta 
much greater than the Fermi momentum, there are no low energy particle or hole excitations. 
Therefore the Green's function should not have an interesting pole structure in the frequency 
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dependence. From the expression for the Green's function (36) we can obtain the leading 
low frequency dependence. Outside of the near horizon Lifshitz region, all quantities are 
analytic as — t- 0. However, the turning point is inside the Lifshitz region. As a; — >■ 0, 



the turning point necessarily has r — )• oo. In this regime we may approximate (27) by 

V = g^ (P - r'^-W) . (42) 

The turning point is therefore 

n = (43) 



Using the above two expressions for the integral in (36) we can easily extract the singular 
frequency dependence of the Green's function. 

Gf(..,A;)«iexp|-2r/7(^y"| , (44) 

where we are working up to an overall momentum-dependent constant and the constant in 
the exponent is 



, = ^^4^. (45) 

2(2-1) 



The expressions (44) and (45) are our final answer for the retarded Green's function in case 
I. The same formulae were derived, for a slightly different question, in [15j. The spectral 
density has no features and is exponentially small at low frequencies but has a branch cut at 
u; = due to the possibility of dissipating (inefficiently) into the quantum critical Lifshitz 
modes. 

Case II is both similar and distinct from Fermi liquids. In a Fermi liquid there is a unique 
momentum hp close to which there are low energy particle and hole excitations. This leads 
to a pole in the Green's function at = and k = kp. In the non- Fermi liquids dual to 
electron stars, the Fermi surface is 'smeared' [21] into a near-continuum of Fermi surfaces. 
These correspond to bulk Fermi surfaces at differing radii in the electron star. This is a large 
N holographic manifestation of the smearing of spectral density due to strong interactions 
with critical bosonic modes. As we might have anticipated, we see that the Fermi liquid 
pole at a single momentum and energy is replaced by a sequence of poles closely spaced 



in energy for any given momentum (40). This is because any given momentum (not too 
far) below k^p is the Fermi momentum at two radii in the bulk. There are therefore gapless 
excitations at these radii as well as arbitrarily light excitations at nearby radii. 

Here we should emphasize that the above results show (unsurprisingly) that a WKB 
treatment is powerful enough to identify a large number of closely spaced yet discrete poles, 
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rather than a continuum. We will have to do some work below in order to re-establish 
contact with the picture of a smeared Fermi surface. 

As we wish only to capture nonanalytic dependences on frequency and momentum, we 



can zoom in on the poles (40) of the retarded Green's function (124). Thus we can write 



the prefactor ( 38 ) as 



vrn 



(46) 



Recall that Y includes a term plus ivr in (39). The real part of the denominator vanishes 



at momenta k^p^ where 



Y{0,kP) 



vrn . 



(47) 



These are the many Fermi momenta. Recall that Y is large in the WKB limit and therefore 
these poles are reliable for large values of n. The poles are therefore closely spaced in 
momentum. We can now expand X and Y and the overall integral in (37) about uj = 



and k = k 



(n) 

p to obtain an expression for the Green's function that captures all of the low 



frequency singular behaviour: 



G 



uyuj,k) oc 



E 



0<4"'<fcF ^ + Vn{k - kP) + icn'y iexp<i-27 



■q (^kP^UJ 



(48) 

Here an,bn,Cn,Vn are order one (in 7) positive numbers that can be obtained from the 
integrals appearing in (37) at = and k = kp . In particular 



d^Y{0,kP) 



(49) 



and 



ri 



\/vio,kP] 



m 
r 



dr + rh log ri . (50) 
Here we added and subtracted a term rh/r'm. the integrand so that we could take the integral 



all the way to the boundary at r = 0. The constant in (48) is again given by (45). The 
answer agrees precisely with the 'semi-holographic' result with a Lifshitz IR geometry, 
as it must. The singular frequency dependence of the imaginary part of the denominator 
follows from the same argument that was outlined for case I above. The numerator has 



no singular frequency dependence because the turning point for the overall integral in (37), 
namely ri, is in the far region where the zero frequency limit may be taken smoothly. 



Equation (48) is our final expression for the retarded Green's function in case II. It de- 
scribes a large number of Fermi surfaces making up the star, with Fermi momenta bounded 
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between zero and the extremal Fermi momentum kp of (23). These fermions are exponen- 
tially long lived at low frequencies. The spectral density will therefore be well approximated 



by a sum of delta functions. Note however that the residue of the poles in (48) depends 
exponentially on 7 times a function of the Fermi momenta. 

The plots in figure [5] show the exponent — a„ of the Fermi surface spectral weight, 



-270,1 



(n) 

in (48), as a function of the Fermi momenta kp for several values of the electron 



star parameters. We see that the Fermi surfaces at small momenta, near the top of the 
potential, have exponentially larger spectral weight than those near kp, at the bottom of 
the potential. We can understand this intuitively from the plot of the effective Schrodinger 
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Figure 5: Exponent of the residue of the Fermi surface poles as a function of Fermi mo- 
mentum (in units of the charge density). Left hand plot has m = 0.3 for all curves while 
the right hand plot has z = 2 for all curves. The plots show that Fermi surfaces near the 
maximal Fermi momentum have exponentially small spectral weight compared the other 
poles. 



potential in figure [3j the fermion has further to tunnel from the boundary through to the 
bottom of the potential compared to high up in the potential. On the other hand, states 

(n) 

near the bottom of the potential are exponentially more stable than those at smaller kp 



according to ( 48 ) . This is because they have further to tunnel to reach the interior Lifshitz 
geometry. Low energy observables can depend on both the lifetime and the spectral weight 
of the Fermi surface excitations. Depending on which exponential dominates for a given 
observable (i.e. lifetime or weight), the dominant fermion contribution will come from very 
different excitations of the electron star. 

Case III is distinct from Fermi liquids. At sufficiently small momenta even the 'smeared' 
Fermi surface is ceasing to exist. The appearance of the IR quantum critical dispersion 
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relation oj ^ suggests an interpretation for the physics at these momenta. Namely, 
when uj > the fermionic excitations are able to dissipate efficiently into on shell quantum 



critical modes and therefore lose their quasiparticle residue. In this case, from equation (36) 
the Green's function has only analytic u dependence. This follows from the fact that ri is 
not in the near horizon region, i.e. that, unlike case I, the fermions do not tunnel directly 
into the Lifshitz sector, and therefore the integrand may be expanded in analytic powers of 
small oj over the whole integration range. Our final expression for the low frequency Green's 
function in this regime is simply 

Gf^^(a;,A:) ocie-2T^ (51) 

Here we suppress analytic momentum and frequency dependence, to leading order the ex- 
ponent is simply a constant c. We should recall however that the boundary of case III is 



itself defined via a non-analytic relation (28) between frequency and momentum: uj > k^ . 

The above analysis breaks down at frequencies and momenta for which turning points 
appear or disappear. These are the boundaries between cases I and II and between cases 
II and III. We will shortly see how interesting frequency dependences emerge at these 
boundaries. Before turning to these points however, we can note an interesting way in 
which the Luttinger count works out from the singularity structure of the fermion Green's 
function. 



5 The Luttinger count 

The Luttinger theorem relates the total charge density of a Fermi liquid to the sum of 
the volumes of Fermi surfaces [30J. For a spin half fermion in 2-1-1 dimensions with Fermi 
surfaces with momentum space volumes Vn 

n ^ ' 

This 'Luttinger count', in which the UV charge density is recovered in terms of the IR Fermi 
surfaces can be violated if the microscopic fermions in the system are coupled to gauge fields. 
A recent discussion of this fact can be found in |31j . As well as Fermi liquid- like states, in 
which all the charge is visible in gauge-invariant Fermi surfaces, it is also possible to have 
'Fractionalized Fermi liquids' [32] in which the sum of gauge-invariant Fermi volumes is less 
than the total charge. It is also emphasized in [31] that the best established examples of 
gauge-gravity duality have the structure necessary to admit Fermi-liquid like (i.e. satisfying 
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the Luttinger theorem) as weh as fractionahsed Fermi hquid phases. Now that we have 
characterised the singularity structure of the fermion correlators of the fermions making up 
the star, we can ask whether the Luttinger count holds true. We can thereby determine 
what type of phase we are in. 

In [23], an argument was given indicating that the Luttinger count did work out for 
electron starsj^ All of the charge is carried by fermions in the bulk. These fermions are 
locally Fermi liquids and therefore at each radius satisfy a Luttinger theorem in the bulk. 
Thus the total charge will be given by the sum over bulk Fermi volumes. This argument is 
not entirely satisfactory, however, because the local bulk Fermi volumes are three dimen- 
sional while the boundary Fermi surfaces are two dimensional. With the boundary field 
theory fermion Green's functions at hand we are now in a position to make a more precise 
statement. 

As shown in [¥1], the charge density of the field theory is given by integrating over all 
the fermions in the star 

Q=^ drr^kl{r). (53) 



^ I s 

The local Fermi momentum kp{r) was defined in ( |21[ ). Note that there are no hats over 
any of these quantities. We have also used the fact that f3 defined in [21] is equal to vr"^ for 
spin half fermions. This expression is compatible with a bulk (three spatial dimensional) 
Luttinger count at each point because we could write it as 

Q = drT^g—^V{r), (54) 

where here V{r) = ^kp{r) is the Fermi surface volume at the radius r. 
From a 2+1 dimensional field theory perspective we need to evaluate 

2 



-.Vr 



(2^)2 ' " 27r 

(n) 

where kp are the Fermi momenta appearing in our Green's function. They are defined 



by poles of the Green's function at a; = 0. From (39) and (47), and more generally the 
discussion of the previous section, we have 

rr2 



f 



dr r^^J kl{r) - k^ ' = n {n + ^) . (56) 
There is no overall 7 on the right hand side because there are no hats on the momenta in 



the square root, cf. (12). The order one contribution to the Luttinger sum comes from 



^This statement should not be confused with an additional statement made in 24 , which is that the 
Fermi surface area as extracted by quantum oscillations would not satisfy a Luttinger count. 
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highly excited states in our WKB hmit and so we can approximate the sum over n in (55) 
by an integral. Letting /c^^^ = E{n) we have 

2 „ 1 /■ , 1 f^P^ ,„„dn 



(2vr)2 7o 



1 r''*F r2 r /g 

I dEE / dr = (58) 

Jri Jkl{r)-E 



1 /-oo (■kpK'T) E 



drr^ / dE = (59) 



1 

3^ 



oo 



drr^tPir) (60) 



g . (61) 



In the first line, kp is the maximum of kp(r) as defined in (23). To get from the first to 
the second line we differentiated the expression ( |56[ ) with respect to n. From the second to 
the third line we exchanged the orders of integration, being careful with limits. We then 



perform the integral over E in the third line. The final result then follows from (53). 

The computation we have just performed confirms that the Luttinger count holds for 
electron stars: the sum of the areas of the gauge invariant field theory Fermi surfaces adds 
up to the total charge density. There is a something slightly magical almost about the way 



that the square root appearing in the standard WKB formula ( 56 ) acts as a projection from 
the spherical bulk Fermi surfaces to the circular boundary theory Fermi surfaces in such a 
way that the factors of vr work out. 

We end this section by emphasizing that while the Luttinger count works out, the dual 
state of matter described by the electron star is not a Fermi liquid. In particular, both 
the spacetime geometry and the fermion Green's function indicate that there is a critical 
Lifshitz sector with gapless degrees of freedom interacting with the fermionic excitations at 
arbitrarily low energies and momenta. 



6 The boundary fermion Green's function: k kp 

This section computes the fermion Green's function at low frequencies and at momenta on 
the boundary between cases I and II above. Recall that this corresponds to momenta close 
to the extremal Fermi momentum kp. This is the largest of many Fermi momenta of the 
electron star. It was shown in [23] that this momentum uniquely determined the period 
of de Hass-van Alphen quantum oscillations of the system. We expect that the Fermi 
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momentum poles will coalesce into a branch cut in the smearing limit. We may expect that 
one endpoint of the cut will be precisely at k ^ kp and therefore we need to characterise 
these poles carefully. 

The WKB approximation used in the previous section fails for momenta sufficiently 
close to kp. The reason for this is that once the turning point becomes sufficiently close to 
the minimum of the potential one cannot treat the potential as locally linear. Specifically, 
near the minimum we can write the potential as 

h2 

V = -A + —{r-rof . (62) 

The radius ro is slightly shifted from as we show in appendix [B| Our discussion will be 
more transparent if we rescale variables and introduce 



y = VjB{r -ro), e = -fA/B . (63) 

The Schrodinger problem then becomes 

- + (^-e + $ = . (64) 

Here dots denote derivatives with respect to y. Let us linearise near the turning point 
and determine whether WKB is valid. Set y = + 6y. To expand and obtain a linear 
potential we require 6y <^ ^/e. Expanding, the equation becomes — + ^/Eby^ = 0. Now 



plugging into the conditions for WKB to hold ( 30 ) , and recalling to convert the derivatives 
with respect to r into derivatives with respect to y, we see that WKB requires e~^^^ <^ 5y. 
From the two inequalities we have just obtained, it follows that we can only consistently 
perform a matching to WKB across a linear matching point close to a minimum of the 
potential when 

1 <C e . (for linear matching) . (65) 

This condition will not hold for k sufficiently close to k*p where e tends to zero. 



When (65) does not hold we will instead solve the Schrodinger equation in the full 



quadratic region where (64) holds and match this solution directly onto WKB on each side. 
The matching solution is now given by Parabolic Cylinder D functions (this is just a fancy 
name for a combination of Confluent Hypergeometric functions) , and not the Airy functions 
of the simpler linear matching. In order for this matching to be valid, we need to be able to 



neglect third and higher order corrects to the potential (62). Generically, recall the rescaling 



(63), this will require y <^ -y/7. Furthermore, in order to be able to match, we need the 
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quadratic form ( 64 ) of the potential to hold simultaneously with the WKB conditions ( 30 ) . 
A sufficient condition for satisfying the WKB conditions together with y <^ ^Jj is that 

e <C 7 . (for quadratic matching) . (66) 



Fortunately we see that (65) and (66) overlap in our large 7 limit and therefore we can 
cover the whole range of momenta by using one or the other method. 



To implement the quadratic matching we need to expand the potential (20) about its 



minimum. At a; = and k = kp, see equations (23) and (24) above, the minimum is at 
r = and V{r^) = 0. We are interested in the potential at nearby values of frequency and 
momentum, uj = 5Cj and k = k*p + 6k. In appendix |b| we show that the potential takes the 



form (62) with 
A = 2g{r^ 



-SCj - k*prl5k 



and 



B = 2nJ -g{n)k*pk'J,{r^ 



(67) 



Recall that r^, is defined via (23) and kp{r) is given by (21 ). From here we can obtain e via 



(63). Before presenting the full Green's function, note that the insofar as we are interested 
in the low lying normal modes and neglect exponentially small tunneling probabilities, we 
can simply assume that the wavefunction is normalisable in the quadratic region. This 
corresponds to the familiar Hermite functions and hence we obtain a spectrum of fermionic 
excitations 



n + 



(n > 0) 



(68) 



For 1 <^ n these are the poles we obtained in our previous analysis (46 ). We have just shown 
that the approach of this current section allows us to go to lower poles n ^ 7. At larger 
values of n the states probe the deviation of the potential away from the purely quadratic 
form near the minimum. 



If we set 6uj = we can think of equation ( 68 ) as determining the Fermi momenta of 



each species of particle in our star. Using the above few equations: 



:(«) 



gir^)rfk*p 



n 1 

7 27 



(69) 



The +^ term here indicates that the extremal Fermi momentum k*p receives a quantum 
correction. As explained below equation (18) above, there are additional corrections at this 
order that we are not considering, so the quantum shift cannot be considered too literally. 
In the large 7 limit the spacing between states is suppressed, leading to a pseudo continuum 



distribution of poles. These poles merge into the poles of (46), with the two descriptions 
overlapping for 1 <C n <C 7. We shall verify the agreement shortly. 
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Similarly if we set 6k = we obtain a spectrum of frequencies Jo) cx ^ + 2^- These poles 
at positive frequencies corresponding to the cost of creating particles of momentum kp for 
the different species of particles at different radii in the star. At kp all Fermi seas of all 
species are empty, so it is not possible to destroy a particle of this momentum. That is why 
there are no negative frequency poles. 

Our discussion so far has neglected the exponentially small imaginary part of the poles. 
We will now compute the full Green's function for k ^ k*p, including the real and imaginary 
part of the poles and their residues. We have relegated the details to appendix [B] Here we 
simply quote the result for the retarded Green's function 

1 _ l^^g-27/s-2eTs 



Recall e was introduced in (63) in terms of A and B which are given in (67) and which 



characterise the minimum of the potential according to ( 62 ) . The overall exponential terms 
characterising the tunneling into the spacetime from the boundary are 

Ib = ^\/Ay- dr + m log To, (71) 

^« = r(-2;lf+F^)*-'°^^''°- '''' 

Here we set V = —A + AV. The first line includes terms that regulate the integral near 
the boundary r — )• 0. The exponential terms that measure the tunneling through the star 
to the horizon, giving an imaginary part to our poles, have exponents 

Is = r^VAVdr, (73) 

J To 

Here recall that is the (linear in this case) turning point closest to the horizon. 



If we ignore the exponentially small imaginary parts in the Green's function ( 70 ) it is 
clear that we recover the poles at e = ^ + n described above. Furthermore, we can 'smear' 
over the poles by looking at the Green's function in the regime 1 ^ e ^ 7. This is the 
regime of overlap of the quadratic and linear WKB matching, so we did not have to do a 
quadratic matching to access this regime. Here the poles should merge and form a branch 



cut. Using Stirling's expansion for the gamma functions in (70) and noting that B is order 
one in 7, one obtains the leading order result 

Ci^bdy.l^,^) ~ ie-2^^B+.iog(7/.)^ (1«£«7). (75) 
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Here # and do not depend on the WKB parameter, the frequency or momentum. Thus 
we see that smearing over the poles does lead to a branch cut emanating from the extremal 
Fermi surface at 

6j ~ v^ik*F - k) . (77) 

The singularity at this extremal momentum is weak. The spectral weight in the branch cut 
increases exponentially as we move to lower momentum fc, away from the branch point. This 
is consistent, as it had to be, with our discussion around figure ^ above; the residue of the 
Fermi surface poles close to k*p are exponentially small. Nonetheless, the extremal Fermi 



momentum does produce a nonanalyticity in the smeared Green's function (76). This can 
then lead to important signatures such as quantum oscillations [24j that pick out momentum 
space nonalyticities and are independent of the spectral weight. The disjunction between 
quantum oscillations and spectral weight is an interesting feature of electron stars. 

7 The boundary fermion Green's function: k ~ oj^/^ 

This section computes the fermion Green's function at the boundary of cases II and III 



above. Recall from equation (28) that this corresponds to small frequencies and momenta 
close to a dispersion relation k ~ w^/^. For momenta below the dispersion relation, the 
fermion states become highly unstable to dissipation into the quantum critical sector. In 
the Schrodinger equation this translates into the fact that we will be considering states 
close to the top of the effective potential. In contrast, the previous section considered the 
most stable states, at the bottom of the potential. Our computation is similar to that of 
the previous section, except that now we must match the WKB region onto an inverted 
harmonic oscillator potential. A further potential complication is that while we are taking 
7 large, as usual, we also need to take a; and k very small as they scale towards the critical 
dispersion relation. 

We will be expanding around the top of the barrier, of for instance figure[3j that prevents 
fermion states from falling into the black hole. Near the maximum the potential takes the 
form 

V = -A - —{r-rof . (78) 



The difference from ( 62 ) of the previous section is the opposite sign of the quadratic term. 



Of course A,B,ro take different values as we shall see shortly. Again rescaling as in (63) 



24 



we obtain the Schrodinger problem near the maximum 

- + (^-e - $ = . (79) 

Also as previously, linear matching close to the maximum of the potential is legitimate when 
1 <^ e while quadratic matching, the topic of this section, requires e <^ j. 

We noted above that this maximum necessarily occurs in the near horizon Lifshitz region 



of the spacetime. The potential we are expanding about the maximum is therefore (27). 
The maximum is seen to occur at 

xo = r^o^ = ^^^^ [h^{2-z)± V/i^z2-4m2(z-l)) . (80) 

where we should pick the plus sign to find a solution for positive Cj. We saw above that 
case III never arises at negative w. We have introduced Xq here; it is simply a number with 
no oj dependence. Expanding the potential about the maximum we can read off A and B 



in (jTSj): 

A = ,00 (^^^^^^^^^^^/^-P), (81) 
B = u^/^^^^^/h^iz - 2)xo + 2(z - l)xl . (82) 

Xo 

Setting A = gives us the precise dispersion relation separating cases II and III. The 
condition for the validity of matching WKB across the quadratic potential is A/B <C 1, 
which is now seen to require that IP' jQp'l^ be close to the critical dispersion relation 

In appendix[C]we perform the matching across the quadratic region to obtain the Green's 
function. The result is 



<-^II-III bdy. Of ^6 ■ ■ • \po) 



Here Ib again gives the amplitude to tunnel in from the boundary 



ri 



m 



Ib = I y VV ) dr + Thlogn , (84) 



^We should additionally worry whether co and k small might invalidate the WKB analysis. The poten- 
tially dangerous assumption is that the quadratic term in the potential dominates all higher order terms 
consistently with the WKB conditions (30 1. We can check that 5r ro = {xo/Cj)^ is sufficient for the 
quadratic term to dominate, while validity of WKB requires \/"^B Sr ^ 1. Compatibility of these two in- 
equalities implies ^J^B{xo/Cj)^ 3> 1. This condition is always true at large 7, the factors of u) drop out. We 
will consider below the ceise when Xo becomes small. 
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while the integrals characterising the periods of oscillation within the well are 



Iw 
Tw 



V-AVdr, 
B 



dr + log y/^B\ro — ri\ 



Recall from ( 63 ) that 



e = 'jA/B oc 7 



(hoc + Xq)^ - 



■,2/z 



(85) 
(86) 

(87) 



The previous section, in e.g. (69) and below, described how the exponentially stable 
fermionic excitations ended at the bottom of the potential in a harmonic oscillator fashion. 



We can now use the Green's function (83) to see what happens to the poles near the 
maximum of the potential. In appendix [C] we find the poles of the Green's function in the 



regime 1 ^ |e| ^ 7. We also derive simple expressions for and Tw in (167) that we use 
in the remainder. In this regime the Green's function ( |83| ) becomes 



III bdy. 



OC 



tan 




log 4 

e uj 



(K e <C 7) • 
(88) 



Here we used the expressions in (167) of appendix [C] for and Iw- Writing e in terms of 
real and imaginary parts 

e = £1 + ie2 , (89) 



the poles are at 



£1 



£2 



27V5^\//i^ - 1 

log - 

z log 7 LO 



7r(2n+ 1) 
log 7 







log 7 



if £1 < 
if ei > 



(90) 
(91) 



We see that, in this regime, while changing the value of (real) Cj shifts the poles around, 
they all lie in a line that tilts at the origin of the e complex plane. Ignoring the explicit Co 



dependence in (90), then setting e = jA/B and solving for 6j maps the line of poles with 
a finite imaginary part into a semicircle in the complex cu^/^ plane. We see this semicircle 
vividly in the plots of section [8] below. The explicit logarithm of oj only leads to mild 



violations of the semicircle form. A sketch of the poles described by (91) in the complex 



frequency plane is shown in figure [6j The plot is illustrative only, for (91) to be valid one 
should take a larger value of 7, which would see the second set of poles tilt only slightly 
into the lower half plane. 
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Figure 6: Typical distribution of the poles in the fermion Green's function close to k 



l/z 



as described by (91). This plot has 7 = 100 and z = 2 and order one values for k,rh, etc. 



We have not worried about carefully normalising the uj axes. 



Unlike in the previous section, we have just seen that the limit 1 ^ |e| ^ 7 does not 
smear the Fermi surface poles into a branch cut, as the poles are still present and given by 



(90) and (91) and plotted in figure 6l Me must therefore work a little more to quantify how 



the Green's function changes at e = 0. The tangent appearing in the Green's function (88) 
can be written as a sum over its poles 

00 ^ 

tanF(e) = - J2 j,, , ^ , ^1, • (92) 
Now let us sit on (or parametrically close to) the real e axis and ask whether we can 



approximate this sum by an integral. For e < the function F{e) in (92) is real to the 



order we are working and the tangent has poles on the real axis at (90). The spectral density 
(imaginary part of the Green's function) is thus a sum of delta functions and we should not 
replace the sum by an integral. Note that while the residue of the poles decreases as one 



moves away from e = 0, due to the overall e"^"^^^ factor in (88), it does not vary strongly 
over the spacing between poles and therefore one cannot keep only the first pole. For e > 
however, F{e) has an imaginary part of order e. We can therefore turn the sum over n into 
an integral without coming close to any of the poles. Performing the resulting integral over 
n by contour integration then gives the constant Green's function oc ie~'^"'^^ for e > 0. 
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It is clear from (87) that e changes sign at the dispersion relation 



Xq 

So thus we find that for Cj > uj(k) the Green's function oc ie~'^'^^^ , in agreement with our 



previously obtained (51). Above the dispersion relation (93) the fermion spectral density 
is featureless. Below it is characterised by closely spaced poles. The Green's function (88) 
therefore interpolates between cases II and III as we should expect. We will also see this 
transition from poles to a constant explicitly in the following 'small star limit' section. 

Contribution to the optical conductivity 

While we are focussing on the fermion Green's function, it is of interest here to point out a 
consequence of the existence of the scaling regime w ~ /c^ for a more complicated observable: 
the optical (i.e. frequency dependent) conductivity. The contribution of the stable fermions 
to the optical conductivity was computed in |21j . It was found that (7(0) ~ Q,^ for small Q. 
Here we will estimate the one loop contribution of the states around the critical dispersion 
relation to see if they can dominate the tree level result of [21] at small frequencies. 

The strategy will be to isolate the critical contribution to the conductivity. Since the 



Green's function, given by (83), has (almost) a scaling behavior, we will be able to esti- 
mate easily the contribution of this regime to the conductivity. We will follow closely the 
discussion in [33], where one can also estimate the conductivity from scaling arguments]^ 
The starting point is the expression for the conductivity at one loop: 

(^i^)^ ^ I dujidu2 ■^(^^)~{(^2) ^2 f dkklmG^(uji,k)lmG^(u2,k) , (94) 
lil J wi — a;2 — 'i — J 

where f{i^) is the Fermi-Dirac distribution at zero temperature given by 1 at w < and 
for a; > 0. A is a radial integral over the interaction vertex and the fermion and gauge field 
wavefunctions [33j. This integral extends, in principle, all the way from the boundary to 
the horizon of the holographic radial coordinate and it depends on the frequencies. Since we 
are only concerned with extracting an interesting contribution for a(Q) we can cut off the 
integral away from the horizon and consider a finite UV contribution where A is analytic 
in Q. Furthermore, in the UV, it is clear that the wave functions have a smooth limit as O 
goes to zero and we can replace A by a constant Ajjv^ ^is in |33j . 



''One important difTerence between the present case and the one discussed in [33] is that, because we have 
Fermi surfaces at arbitrarily small momentum k^p\ leading to the scaling form (83 1, we can scale k with O, 



while in [33| the k integral has to be dealt with before one can make a scaling argument. 
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The Green's function G^{u!,k) in (118) is in fact just a function of e and lyy and T\y. 
In appendix [c] we show that Tyy is a 7 dependent constant, while Iw behaves as loguj^ in 
an appropriate limit. Let us then perform the following rescaling: 



(95) 



Under this scaling e remains unchanged (87). This means that the Green's function is of 
a scaling form G(k'^ /co^), up to log-periodic violations which we will neglect. The function 



f{u>) is also invariant under rescalings. Then, we can write (94) as: 

/(wi) - f{u}2) 



Ol — 0)2 — 1 



le 



A2 



dkklmGl k'^/ul''\ Im G ik'^/Cjl'" 



(96) 



The conclusion is that the contribution of critically scaling fermions to the optical con- 
ductivity dominates the contribution of stable Fermi surfaces in the small frequency limit. 



8 The 'small star' limit 

In this section we consider a limit in which essentially the entire star is contained within 
the Lifshitz region. Because the potential is known analytically in this case, we are able to 
characterise aspects of the Green's function explicitly. This limit is achieved by taking z 
and rh such that the mass is only just below the critical mass above which the electron star 
ceases to exist ([T)). Namely 



A' = /i^ - < 1 . (97) 



We expect that in this limit the electron star will be 'small', in a sense which we will now 
make precise. This limit was considered in )21j . as the large /3 limit 

Let us think of keeping m fixed and expanding in A ^ 1. Recalling from ([6]) that 
/i^ = we have 

1 A^ A^ 

2 = + T\ ^2\2 ' /ioo=m-F — H . (98) 

1 — (1 — ra'^Y 2m 

The remaining constant characterizing the IR Lifshitz geometry is computed from the for- 
mulae in [21] to be 

_6-7m2-K2m4 {6 - Im"^ + 2m'^){l + Arn^) 2 

- 6(l-m2)2 + 12m2(i_^2)3 ^ + ' ' ' ' (^9) 



""The quantities are related by A = . ^ — 7 -2^5)2/3 ^273"- is clear that using A is more convenient 
for our purposes here. 
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Also from [21] we can recall that moving away from the IR Lifshitz region (r — t- oo) is 
characterised by an exponent \a\ 

The exponent was computed in |21j and in the small A limit becomes 



(100) 



_ 3i/2m(2-m2)i/2 i i 

(i_ ^2)1/2 a"-^" 2(l-m2) ■ 

The observation to make here is that the exponent |a| is becoming large as A — t- 0. This leads 
to an important simplification in the structure of the star. It means that the corrections 
to the IR Lifshitz behaviour are exponentially small for r > 1 (we will see shortly that the 
coefficients fi,gi,hi can be taken not to be also exponentially large or small in A). Thus 



for r > 1 the Schrodinger equation (18) takes the following Lifshitz form to leading order 
at small A 

^„ off) A2 + 2mr^a) + r2^a;2 \ 

= f [k' J <^ , (102) 

where for convenience we introduced 

f = I'at^ ■ (103) 

In this last expression corresponds to the leading coefficient in the expansion of Qoq 
above. 

In fact, more is true. The star ends also at r = 1. This can be seen numerically or, 
alternatively, as follows. By solving the equations of motion [21] perturbatively one finds 



that the coefficients appearing in the correction ( 100 ) can be taken to leading order in small 
A to be 

9 , 2-3^/2^3n _ ^2)1/2 ^ 

/, = -l, gi = -2m\ hi = -y (104) 



The star ends at the radius r = Vg where the potential in the Schrodinger equation (18) 



vanishes at Cj = k = Q. Using the first order corrections of (104) it is easy to find that to 
leading order 

-Ti = — • (105) 

M 2m2 1 - rn2 ^ > 



From this result we learn two things. Firstly that the correction in (100) is still small at 
the boundary of the star, and therefore our computation of Vg is selfconsistent. Secondly, 
to leading order we have = 1. Thus we find that the small A limit of electron stars has 
an interesting and simple structure: for r > 1 (inside the star) it is Lifshitz and for r < 1 
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(outside the star) it is Reissner-Nordstrom-AdS. This is potentiahy physically exciting as 
it means that all the fermions in the star are in the quantum critical geometry. 

The local fermion density is therefore discontinuous at the edge of the star in the strict 
small A limit. At r = 1 the density abruptly falls to zero. This kind of sharp endpoint of 
a 'Fermi surface-like' object in the bulk is familiar from \34:\ 135) where instead of fermions 
there are endpoints of strings ending on a D brane and electrostatically repelling each other. 
A similar jump from Lifshitz to Reissner-Nordstrom is discussed in section 7.3 of [20]. The 
reason that these discontinuities can arise is because we are looking for solutions subject to 
the constraint that there is a charge density at the boundary. 

The fact that the boundary of the star occurs at = 1 is not a coordinate-invariant 



statement and ultimately depends on the choice of normalisation we made in ( 104). We will 
later want to plot only physical, dimensionless ratios of quantities in which this ambiguity 
will cancel out. In order to do this, we need to compute the chemical potential of the 
dual field theory, /t, as well as the 'speed of light', c. This can be obtained in the small 
star limit as followsj^ Firstly note that the definitions of fi and c appearing in Q can be 
extended away from the exterior Reissner-Nordstrom- AdS region by defining c(r) = r^\fj9 
and c/i(r) = h — rh' . The boundary values are then c = c(0) = c(rs) and fi = /i(0) = fi{rs)- 
Diff'erentiating these expressions, using the background equations of motion studied in |21j 
and then integrating, one obtains 

/"OO /* OO ^ 

c=\ r^g^/^adr , cfi = "7? (/ " h^h'h) dr . (106) 

In the small A limit all of these integrals may be evaluated on the Lifshitz background and 
with Ts = 1. One obtains 

^=2(^' ^ = ^['^^)- 

Here the charge density a is an order one constant in the small A limit. Using the definition 
of a in [21 J one finds 

3'^m(l — m^) 
^=2rn4-7m2 + 6+---- 

To start with it is instructive to consider the bound state spectrum. These are bound 



states of the Schrodinger problem (102) with cd = 0. Solving this equation will determine 
the Fermi momenta kp of (56) explicitly in the small star limit. At zero frequency it is 



This approach originated in collaboration with Pavel Petrov and was used to prove some assertions made 
in |36l- Similar methods were also used in 1371. 
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easy to solve the Schrodinger equation subject to normalisability in the interior r — t- oo to 
obtain the modified Bessel function 



Here the coefficient 



(109) 



(110) 



The bound states are just given by imposing <l> = at r = 1, as the potential is exponentially 
large (at small A) at r < 1. The zeros of the Bessel function can be found explicitly in the 
small momentum limit 7A: <C 1 to be 



(111) 



We see that the Fermi surfaces are accumulating exponentially at low momentum. This is 
because the Schrodinger potential has a maximum at the origin of the star. For these states 



to exist we need E M and thus, from (110), 27A > 1 . 



Solvable case: z=2 



At finite frequencies the Schrodinger equation ( 102 ) is still not solvable in general. However, 



it is clear that z = 2 is a special value as only three different powers of r (as opposed to 



four) appear in the Schrodinger equation (102). When z = 2 the solution to (102) satisfying 
ingoing boundary conditions at the horizon is the Tricomi confluent hypergeometric function 



+ 



- 2mcu)7 



(112) 



Where v is again as in (110) and real. When 2; = 2, to leading order at small A, the mass 



m 



1/2, from e.g. (98). The quasinormal modes are found by imposing ^*(1) = 0. This 



condition can be solved explicitly in the low frequency limit 70) <C 1 to give 



(-Z7w)" 



l+iu _|_ i-^{k'^-2mCj) \ 



46j 



r ( + r (ii.) 



(113) 



The first computation we can perform with this result is to derive the dispersion of the 
fermion modes about the Fermi momenta (111). Taking the limit 7|(fc^ — 2mcd)/u\ ^ 1 of 



the above formula, one obtains the dispersions 



k'-k 



{n)2 



2mu) = . 



(114) 
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This result is self consistent so long as 7^A:^''^ ^ ll'^l- We have therefore found the Fermi 
velocities of these poles. 

The next step is to find the imaginary part of the dispersions. In fact simply by expand- 



ing (113) it turns out that the pole (114) never acquires an imaginary part at any order in 
the expansion in 7|(/c^ — 2muj)/Ld\ ^ 1. This remains true if we keep higher order terms in 
the expansion of the confluent hypergeometric function in juj. To extract the exponentially 
small imaginary part we need to be more careful with the asymptotic series expansions of 
the gamma functions. In particular, we need to keep nonperturbative corrections to the 
usual expansion about Stirling's formula. Following [3H] we can note that if we write 



T{z) = V2^z'-h-'e^^'\ (115) 

then the reflection formula for gamma functions implies that 

,(.) = ?<£>^-l,„g(l (116) 

The first, antisymmetrized, term on the right hand side is the standard asymptotic expan- 
sion in terms of Bernoulli numbers. The second term gives corrections that are exponentially 



small in the upper half plane. Using the above formulae to refine the expansion of (113) we 
find that the dispersion including the leading imaginary part is 

P - kP'- 2rhco + ,4")2!i^e--^^^"^^/2|^l =0. (117) 

Thus for small electron stars we obtain explicitly the exponentially long lifetime that we 



discussed more generally around (48) above 



As well as exponentially stable poles, the expression ( 113[ ) captures the poles that are 



away from the real axis, some of which we discussed in section[7] above. A virtue of the exact 



solution (112) is that it interpolates between cases II and III without needing to perform 



separate WKB matching procedures. In particular, it captures all the poles in the Green's 



function in the single expression (113). This expression is similar to the poles of formula (83) 



for the boundary between cases II and III, except that (113) has a wider regime of validity. 
To make these statements precise, we proceed to find the boundary Green's function. 

To find the Green's function, we need to match the solution of the Dirac equation near 
the boundary of the star onto a solution in the exterior region. To connect with the rest of 
this paper, we should do this in the strict WKB limit. As we discuss in appendix [D] this 



involves taking ~ 7A ^ 1. Expanding (112) as previously in rju <^ 1, using (119) and 



(120) to match the resulting solution onto an exterior WKB solution and then using (35) 
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to read off tlie Green's function, we obtain: 



Am 



y( l + if I ij(k-'-2rh,uj) \-p( ■ s 



(118) 



l — iu _|_ i'^{k^ —2rhCj) 



T{iv) 



We see that the poles of the Green's function essentially agree with (113). The extra factor 



of e*'^/^ is the usual shift in the zero point energy that, while subleading, is captured correctly 
in a WKB treatment. The boundary condition we used above, setting the wavefunction to 
zero at r = 1, was too crude to capture this correctionj^ Perhaps curiously, the Green's 



function (118) bears some resemblance to the IR Green's functions of extremal Reissner- 



Nordstrom-AdS [2]. In the prefactor of (118) we used the fact that uj^k <^ jl in the regime 
we will be interested in. The exponent in the overall suppression e~'^'^^^ due to tunneling 
from the boundary is again given by (84) with ri = 1. The frequency and momentum 



dependence of this term will be weak compared to the remainder of the Green's function 
and so this term is taken to be a constant in the plots below. 
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Figure 7: Poles in the retarded Green's function: Imaginary part of the retarded Green's 
function, lm.G^{Cj, k), as a function of complex frequency at fixed real momentum. Com- 



puted using (118) with rh = l/-v/2) corresponding to z = 2 as required. Left plot has 



7 = 1000, = 40 and k/ji = 2/3000. Right plot has 7 = 10, = 4 and k/fi = 2/30. 



In plotting the Green's function we should use dimensionless variables. For z = 2 we 



find from the formulae in the previous subsection, such as (107), that c = and fi = 3/2. 

'^A linear WKB matching across this turning point means we cannot access near extermal Fermi momenta 
between cases I and II. However we still capture the crossover from cases II to III, our primary interest. 



34 



Plotting the Green's function against (b/{cjl) and k/jl will undo the coordinate choice we 
made in taking = 1 while also fixing the normalisation of time of the boundary quantum 
field theory so that the speed of light is unity. The result is shown in figure [7] 

In figure [7] we see several features discussed previously. There are poles along the 
real frequency axis corresponding to excitations of the many different Fermi surfaces, but 
beyond a) ~ A;^ the spectral function is constant. The poles away from the real axis form 
the semicircle anticipated in section [7j To understand the semicircle we can relate the 



poles of the small star Green's function ( |118 ) to the general Green's function (83) close 
to the dispersion relation w ~ /c^. To connect the two we must take the large v limit 



of (118) keeping e fixed. Of the two more complicated gamma functions appearing in 
the denominator of (118), the one in the numerator is r(^ — e + iu) while the one in 



the denominator is r(i — e). At large v with e fixed we land on the general result (83). 



Repeating the arguments of section [7] gives a semicircle of poles in the lower half plane. 

The semicircles form essentially a right angle with the real axis in figure [7j different from 
the small angle of section [7| In the small star limit the angle is determined by log 7 A ~ log v 



rather than the log 7 of equation (91), see appendix Dl For the values of the parameters 



chosen in the plots, this logarithm is not sufficiently large to lead to a small angle. 



Because the Green's function ( 118 ) captures both cases II and III, we can plot ImG^(c2;,A;) 
to produce a precise version of the sketch of figure |4j The result is shown in figure [Sj 
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Figure 8: Spectral density of many Fermi surfaces interacting with a dissipative critical 



sector: lm.G^{Cj, k), as a function of frequency and momentum. Computed using (118) 
with rh = l/\/2, corresponding to z = 2 as required. The plot has 7 = 100, v = 30. 
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As a final plot, we illustrate more explicitly how the spectral weight of the stable Fermi 



surface poles are dissipated away as they pass through the critical dispersion relation ( 93 ) . 



Figure [9] shows the spectral density of a pole along the real frequency axis as the momentum 
of the pole moves into the unstable region. 
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Figure 9: Dissolution of the Fermi surface pole in the critical region: Spectral density as a 
function of frequency for various values of momentum, from k/jl = 0.053 for the leftmost 
peak to k/fi = 0.08 for the rightmost. The plot has 7 = 10, z/ = 3. 



9 Discussion 

In this paper we have characterised the Green's function G^{uj, k) of gauge-invariant fermion 
operators in phases of 2+1 dimensional field theories at finite charge density dual to electron 
stars. Three main results are 

• G^(0, k) shows a large number of closely spaced Fermi surface for k < kp. The 
volumes of these Fermi surfaces add up to the charge density, thereby satisfying the 
2+1 dimensional Luttinger count. 

• G^{u},k) changes qualitatively at a specific dispersion: uj ~ k^. For real lo < k^ the 
Fermi surface excitations are very stable, while for uj > k'^ they are unstable and the 
spectral density eventually becomes featureless. 

• Smaller kp Fermi surfaces have exponentially larger spectral weight than the more 
stable large kp Fermi surfaces. 

We also estimated that the critical fermions near the dispersion oj ^ k^ may dominate some 
physical observables such as the optical conductivity. 
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As emphasised in [251 [261 [23] , electron (and neutron) stars in the bulk should be thought 
of as Fermi surfaces of generalised free fields. Generalised free fields (see e.g. [39]) are 
the mechanism whereby holography is able to give a semiclassical description of strongly 
interacting physics: the spectral weight of an operator in the field theory is distributed 
among an infinite tower of generalised free fields, which correspond to different harmonics 
of the radial quantm number in the bulk. Thus, the large number of closely spaced Fermi 
surfaces can be thought of as a 'smeared' Fermi surface of the strongly interacting fermion 
operator |24j. We have seen explicitly in this paper that WKB resolves the individual Fermi 
surfaces. Whether or not the spectral weight appears smeared or discrete will depend on 
the observable. 

The fact that the field theory Luttinger count works out for electron stars indicates 
that they are not dual to 'fractionalised' Fermi liquids in the sense of |32l I31|. Nonetheless, 
there is no indication that the field theory SU{N) gauge symmetry is Higgsed in these 
backgrounds, implying that the transverse gauge field excitations should remain gapless 
at the lowest energy scales. The IR Lifshitz geometry of the electron star should then be 
thought of as a gauge-invariant and strong coupling cousin of Landau damping of the gauge 
field excitations by the density of fermions. Many of the excitations of the Fermi surfaces are 
unable to dissipate efficiently into the critical Lifshitz modes. The appearance of a particular 
dispersion uj ^ k^, separating regimes that can and cannot dissipate, has the flavor of a 
'non-relativistic lightcone', suggesting an emergent kinematic constraint controlling the IR 
interactions between bosons and fermions. We feel that this dispersion relation remains to 
be further understood. 

It is of interest to find gravity duals to fractionalised Fermi liquid phases. These will 
be characterised by the fact that the total volume of bulk gauge-invariant Fermi surfaces 
will add up to less than the total charge density (and that there is at least one such 
Fermi surface). This indeed occurs if one neglects the backreaction of the bulk fermions 
on the geometry and studies fermion correlators in the extremal Reissner-Nordstrom-AdS 
background as in |12| [T3l 114] . However, even for the fermions considered in those papers, 
where the effects of backreaction are subleading in the bulk classical limit, it is plausible 
that once backreaction of the bulk Fermi surface is incorporated then all the charge is 
carried by fermions rather than the black hole horizon [20]. This will then likely lead to 
the Luttinger count being satisfied, as we found in this paper for the electron star. The 
conclusion, in that case, would be that the 'large A^' and low energy limits do not commute 
in terms of diagnosing the phase of matter. The electron star limit was designed to have the 
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virtue of capturing the correct low energy physics at leading order in the large A'^ (i.e. bulk 
classical) expansion |21| . The objective is therefore to find a context where a backreacting 
bulk Fermi surface can coexist with a charge-carrying horizon at zero temperature. The 
recently discussed Fermi surfaces in dilaton gravity [IQ] could be one interesting case to 
consider. 
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A Intermediate matching formulae for cases I, II and III 



To start, recall the well known matching conditions across a generic turning point where V 
vanishes linearly in radius. Then, the Schrodinger equation is solved in the matching region 
by Airy functions. Schematically, with x changing sign across the turning point, we have 



1 _T.3/2 

-e ^ 
2 



sm 
cos 



(-x)3/2 + vr/4) 
-x)3/2 + vr/4 



(119) 
(120) 



As we noted in the discussion below equation (18), from the start we are neglecting terms 



that have effects comparable to the 7r/4 phase shifts here. Nonetheless we will keep track of 



the phase shifts with a view to solving the WKB problem posed by equation (18) precisely. 



At the horizon we must impose ingoing boundary conditions. Therefore, in the region 
between r = oo and the first turning point {r^ in cases I and II, ri in case III) we must 
have when uj > that 



r>ri/3 



exp <^ij I dry 

''"1/3 



-V 



(121) 



Using the connection formulae ( 119 ) and ( 120 ) we find that on the other side of the matching 
point 



r<r 



1/3 



e ""'"/^ I exp < 7 



'■1/3 I I 

drVV r + 2 ^-^P 



r'"i/3 
"7 / drVV 



(122) 
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In the cases I and III this region beyond the first matching point is already the outer region. 



Comparing with expressions (31), (33) and (34) we can conchide that the retarded Green's 



function at leading order in the WKB expansion is 

G^{u}, k)oc- hm r'^""^ exp I -2-f / drVv 

2 r->0 



(cases III and I) . 



(123) 



This was formula (36) in the main text. 



In case II we need to propagate the solution through the intermediate regions. Using the 



connection formulae (|119|) and (|120|) two additional times, it is straightforward to obtain 

(case II) (124) 



G^iio, k)(x— lim r"""^ exp { -2-f / dr^V 

2 r-^O I ./^ 



iQ 



-2mL , 



where we specify Q in the formulae below (37) in the main text. 



B Matching formulae for A; ~ 

This appendix contains the computation of the matching across a minimum of the potential 



necessary for section 6l Firstly, we should derive equation (67) which, combined with equa- 



tion (62), describes the potential (20) near the minimum. To zoom in on the problematic 



region near r = r^,, we Co = dCo , k = k*p + 5k and r = r^, + 5r. Expanding the potential 



( 20 ) in all three variables gives (to first order in the frequency and momentum and to second 
order in 6r) 

V{r) = 2 {kyg{r,)5k - j{u)5u:^ + 2 {kyg'{r,)5k - j'{r,)5Cj) 5r 

+ (k*pg"{n)6k - j"{r,)56j - k*pk'J,{r^)g{r,)) 6r^ + ■ ■ ■ . (125) 



In order to keep the length of the formula down, we defined g{r) = r'^g{r) and j{r) = 
g{r)h{r) / f{r). From this last expression we see that turning on the variations Suj,6k moves 
the extremum from 6r = 0. The new extremum is at 



kyg'{r^)6k - j'{n)6oo 
g{r^)k*pk'l,{r^) 



(126) 



to leading order in 6(jj,6k. We wish to follow the behavior near the extremum, so we will 
consider this an expansion around the actual minimum = r^, + 5r^. Substituting into 



(125) and re-expanding about the true minimum the potential becomes 
Vir) = 2 (kyir^)6k - j{n)6uj - ~g{rSkFkUn){r - r^) 



+ 



(127) 
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Recall that since kp is a maximum of kF{r), kp{r^) is negative. Comparing with the 



quadratic potential ( 62 ) in the main text we can read off the values of A and B quoted in 



(67). 



In the main text we then performed the rescaling (63) in order to put the Schrodinger 



+ ^)$ = 0. (128) 



equation in the form (64) which we rewrite here for ease of reference: 

4 

We can now solve this potential exactly and describe how we use the exact solution to 
perform matching. The general solution is 

Hy) = CiD_.^^{y) + C2D_.^^{iy) , (129) 

where Ds{y) are the parabolic cylinder D functions of order s. The asymptotic expansion 
of these functions can now be used to match onto the WKB solutions. The common 
region where both the WKB and the quadratic potential approximations are valid is when 
max(l, y^) <C \y\ <C ^/7. Expanding for large \y\ we can quote the asymptotic expansions 
of 

$(y),^+oo Cie-Ty-|+-+(e— ie-*fC2)e+'^y-5-% (130) 



2 



+ (f^^^i + e-f e^f C2 j y-i- . (131) 

Now we should match these results with the expansions of the WKB wave functions. 
The WKB wavefunction for 5r = r — ro > 0, that is, towards the star interior, takes the 
form 

(72^)3 (72^)3 

In order to evaluate the integrals in the exponents here, it is useful to write, recalling the 



form of the potential (62), V{r) = —A + AV{r). We can then note that given that Tq + 5r 
is in the overlap region {^/A ^ (5r ^ 1), then we have A <^ AV{r) over the whole range of 
integration. Expanding to first order in A we can then write 

/ I- \ 9 A 

To 

5r 

(133) 

Note that a potential logarithmic divergence cancels between the second and third terms in 
the integral. 



Jr,+5r Jro V 2VAV B\r - ro\J 4 S ^ < 
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Using the rescaling (63) we can then express the WKB solution (132) for 5r > in the 
matching region as 



V2R+ 



V2R. 



-fls-sTs 



ilB) 



y 2 "=6 4 



(134) 



where 



Is = \/Aydr, 



B 



+ 



1 



dr - log V75|r3 - ro\ . 



(135) 
(136) 



2VAy 

The subscript 5 is supposed to remind us these are integrals calculated over the repulsion 
generated by the star. Notice that these quantities are not a function of our Schrodinger 
energy e. 



Equating (130) and (134) 



C2 



V2R^ 



(137) 
(138) 



This is one half of the matching procedure. We can now go through the same exercise on 
the other side of the minimum of the potential. The solutions in the matching regime for 
6r = To — r > are: 



V2L+ 



+ 



L. 



y\ 2+^e 4 + 



where now 



ilB) 

Ib = 

Tb = 



V2L_ 

{lB)l' 



-flB-sTB 



y\ 2 =e 4 



(139) 
(140) 



B 



+ 



dr — log \J'~^B\ro\ ■ 



(141) 
(142) 



2^/AF 

Here the subscript B is reminding us that these are integrals to the boundary of AdS. One 
small comment is that integrals diverge near the boundary; rather than integrating to r = 



we should introduce a UV cutoff as in (31 ). Ultimately, in the Green's function, the integral 



will be regularised as in (37). We incorporated the regularisation terms in the corresponding 
equations in the main text. 



41 



Matching the result (140) with (131) we obtain 



L_ 



V2 
^/2 



Co 



27r 



(143) 
(144) 



FinaUy, we can now eUminate the matching region altogether by using the above expression 



together with (137). The coefficients of the WKB solution on either side of the matching 



region are thereby related as follows: 



r(i + s) 



--lIs-eTs 



L. 



(145) 
(146) 



Recall, using the linear matching conditions (119) and (120), that imposoing ingoing 



boundary conditions at the horizon and matching across the turning point r-^ nearest the 
horizon will set R- = Thus we have all the ingredients to find the Green's function. 



Taking r — )• in ( 139 ) we have 



„2mL 



G^(a;,A:) oclimr' ^ 



(147) 



The regularising r^™"^ can be incorporated into the definition of Ib and and we have 
done this in the main text. From our above results 



1 _ 1 „-2^I^-2eT^ 

2r(|+.)^ 



V2n _ 1 „ienp-2"/Is-2eTs 

r(i-.) 2^ ^ 



(148) 



This is the Green's function we have quoted in equation ( 70 ) in the main text 



C Matching formulae for k ^ u 



l/z 



The Schrodinger equation in the regime across which we wish to match is 

-^)^> = 0. 



(149) 



This is exactly the same problem as we have solved in the previous appendix, equation 



(128), but with y — )• \/ty and e — )• —ie. Therefore, the solutions are 



(150) 
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It is important to be careful about the asymptotic expansions here, as the rotation in the 
complex plane does not allows us to just analytically continue the expansions of the D 
functions we used before. In the present case 

/ e 8 e ^4 v27r 



C2 - e"^ 4 e* 8 Ci e * 4 y 



+ 



J/— > + oo 



4 e 



-^Te^^Ci ) e 



+ e^4e*8C2 



_ IT -TtT 

27re ^4 8 



Ci I e+*'ry-5+*" 



(151) 



(152) 



For positive y we need the wave function to match onto a mode that is falling into the 
horizon. This means the coefficient of the e * 4 term needs to vanish. This requires 

•TT _ TT /—— 

e 86 ^4V27r 



C2 



r(i-..) 



Cie"4e*- 



0. 



(153) 



From this we obtain (up to an inessential normalization) the form of the wave function at 
negative y to be 



*(-y) 



27r 



(154) 



V r(i-i.) 

If we were looking for eigenstates with purely outgoing boundary condition we would obtain 
e = \ + n), the analytic continuation of the standard harmonic oscillator problem. Here, 
however, we need to match this wavefunction into a WKB approximation that captures the 
fact that particles can oscillate inside the potential well. The WKB solution takes the form 



(-7V) 



-Vdr 



+ 



Q- 



(155) 



To match this WKB result onto the expansion (154) we use the analogous representation 

6r 



to that in ([133]), letting V = -A + AV. To first order in A 

rVo—Sr 



n 



A 



A 



B{5rf A 

, , , , or log ■ 

2V-Ay B\r-rQ\J 4 B Vo-n 



where 



Recalling the rescalings (63), matching leads to 

ie^^\/27r 
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27=AF \r - 
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(158) 
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where W is supposed to remind us that we are integrating over the well. 

The final step in computing the Green's function is to match the oscillating WKB 



region onto the asymptotic regime. This is a standard linear WKB matching, using (119) 



and (120). From the definition (35) we obtain 



G'« = AiQ_Je-27/s (160) 

2m 2 0+ + i ' ^ ' 

Q— 

where Ib is again the regularized integral characterizing the amplitude to tunnel in from 
the boundary. It is given by 



ri 



m 



Ib = I yVV ] dr + rhlogri . (161) 



The above Green's function is the result we quote in the main text in ( 83 ) . 



It is instructive to identify the poles in the Green's function (160). These occur at 



r |_2 - «ej 

There are two types of solutions to this equation. Most of the poles can be captured using 
Stirling's approximation for the gamma function, i.e. taking |e| ^ 1. Explicitly separating 
into real and imaginary parts: 

e = El + ie2 , (163) 
and taking the logarithm of equation for the poles, we obtain 

ei (vr - arg(e)) - £2 (log|e| - 1) + i [ei (log |e| - 1) + £2 (vr - arg(e))] 

= i7r(2n + 1) + i2'ylw + i2eTw ■ (164) 

Consider the real and imaginary parts of this equation. Firstly we can see that there are 
solutions with £2 = and £1 < if Iw and Ty/ are real. These poles are therefore along 
the negative real £1 axis and have 

£i (log |£i| - 1 - 2Tw) = T^{2n + 1) + 2-iIw , (165) 

where n is constrained to run over all integers such that £1 < 0. Working a little harder 
with the expansion of the gamma function, and we do so in section [8j one can show that 
these poles in fact acquire exponentially small imaginary parts. These poles are the highly 



stable poles just below the top of the potential barrier and they connect onto the poles (48) 
of case II that we discussed in the main text. To make the connection precise one should 
map from £ to a) using e = ^AjB are the expressions for A and B in the main text. That 
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these poles are exponentially close to the real axis justifies the assumption above that Iw 
and T\Y are real. 



A second set of poles is obtained with ei > 0. These poles are forced by (164) to have 
an imaginary part. For simplicity, let us look for poles with ei ^ £2- The real part of the 



poles are again given by (165) but now there is an imaginary part, from solving the real 



part of equation (164) given by 



£2 



eivr 



(166) 



log|ei| -l-2riy 

We will shortly show that Tw is large in the regime we are interested in (1 <^ \e\ ^7), and 
therefore our assumption ei ^ £2 is consistent. In fact, we will find 



Iw 



/go 



^Jhl^-rn^ log — , ri4/~log/7. 



(167) 



Thus putting the two sets of poles together, the real and imaginary parts are 

7r(2n+ 1) 



£2 



27V5^\/^^ - 1 1 

. log - 

z log 7 a; 

if £1 < 

ifei>0 

log 7 ^ 



log 7 



(168) 
(169) 



The geometry of these poles is discussed in the main text around ( 90 ) . 



Finally, let us derive the expressions (167) for T^y and Iw- Consider Iw first. The 



upper endpoint of the integral expression (158) is inside the near horizon Lifshitz region. 
We can split the integral into a part that is contained in the Lifshtiz region and a part that 
is not. Schematically 

(170) 



The precise radius where we choose to make the split will not be important. We are aiming 
to isolating a contribution to the integral that is singular in the limit cj — )■ 0. Away from 
the Lifshitz region, the potential does not have any singular cD dependence. Focussing on 



the near horizon contribution we can change variables to x = r^Cj and obtain, using (27) 

f9^ 



Iw 



'° dx 

X 



^((/ioo+x)2-m2)(l-(x/Xo)2/-). 



(171) 



Recall that Xo is given by (80) and does not depend on uj. The upper limit of the integral 
is therefore regular. The lower limit leads to a logarithmic divergence with cD that we have 



quoted in (167). The precise radius at which we evaluate the lower limit appears inside the 
logarithm and is therefore subleading as c<j — t- 0. 
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For T\Y we can also perform the same change of variables in (159). The integral part of 
Tw is regular near the upper limit Tq by construction, so the same reasoning applies and we 
don't need to worry about this contribution. The lower limit gives us a finite contribution. 
This contribution is subleading with respect to the extra piece of T\y, outside the integral, 
in the large 7 limit. In this terms we see that the uj contributions cancel and there is no 
frequency dependence 

THz-logv^ko-ril -logv^lxoM^/^-logVT- (172) 



This is the expression quoted in (167). 



D The WKB limit for small stars 

In this appendix we comment briefly on the validity of the WKB limit for small stars. 
The potential issue is that while 7 is being taken large as always, we are also taking the 



parameter A of (97) to be small. 

Firstly recall from footnote [H] that in terms of the quantity /3 of f21] we have A ~ /3~^/^. 
Furthermore, recall from |21j that /3 ~ e'^7^. It follows that 7A ~ 7A ~ 7V3g-2/3 ^ 
because we are taking 7 1 while the electron star solution requires e ^ 1. Thus the regime 
of validity of out computations always has 7A ^ 1. We will see that this is sufficient to 
force us into the WKB regime even in the small star limit (and, in particular, is compatible 
with the WKB hmit). 

The most dangerous place in the potential in the small star limit is the boundary region 
between cases II and III that we considered in section [7] in the main text. Recall from 
footnote [3] that a condition for WKB matching to be possible was that y/'~^B{xo/Cj)^ ^ 1. 
Let us check whether this condition is compatible with the small star limit. Using the 



formulae (82) and (80) for B and Xq respectively, and expanding for A ^ 1, we find that 



-fB{xJu:)-^ } ' . (173) 

[ (7A)i/2 z>2 

The discussion in the previous paragraph showed that 7A ^ 1 remained large even while A 
was taken small. This is then sufficient for the WKB condition to hold. Thus we find that 
the small star limit of electron stars is compatible with our WKB analysis. 

Factors of A <C 1 appear in various physical quantities in the small star limit. The 



critical dispersion relation (93) becomes a) ~ A ^k^ for z > 2. Also for z > 2 all the factors 
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of 7 in the locations of the Green's function poles in (90) and (91) become replaced by 7A. 
This allows one to tune the spacing between the poles independently of the scale set by 7. 
It is interesting to calculate Iw and Tw in the small star limit for z = 2. Here we 



can obtain exact expressions that can be matched against (167) where we extracted the 



dominant contributions. We can use expressions (85) and (86) and plug in the exact WKB 



potential AV at z = 2 and k obeying the (exact) critical dispersion relation cj = :^/j^2a ' 

In the small to limit we can perform the integrals exactly and obtain: 

Iw = \/5^A log J- Tw = log VtA , (175) 



as predicted by (167), including the extra factors of A that we mentioned in the previous 
paragraph. 
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